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Abstract 

Background: Tlie microbiota of tine nares lias been widely studied. However, relatively few studies liave 
investigated the microbiota of the nasal cavity posterior to the nares. This distinct environment has the potential to 
contain a distinct microbiota and play an important role in health. 

Results: We obtained 35,142 high-quality bacterial 16S rRNA-encoding gene sequence reads from the nasal cavity 
and oral cavity (the dorsum of the tongue and the buccal mucosa) of 12 healthy adult humans and deposited these 
data in the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI) (Bioproject: 
PRJNA248297). In our initial analysis, we compared the bacterial communities of the nasal cavity and the oral cavity 
from ten of these subjects. The nasal cavity bacterial communities were dominated by Actinobacteria, Firmicutes, 
and Proteobacteria and were statistically distinct from those on the tongue and buccal mucosa. For example, the 
same Staphylococcaceae operational taxonomic unit (OTU) was present in all of the nasal cavity samples, comprising 
up to 55% of the community, but Staphylococcaceae was comparatively uncommon in the oral cavity. 

Conclusions: There are clear differences between nasal cavity microbiota and oral cavity microbiota in healthy adults. 
This study expands our knowledge of the nasal cavity microbiota and the relationship between the microbiota of the 
nasal and oral cavities. 
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Background 

Bacterial communities play important roles in the health 
of their hosts, including roles in immune system devel- 
opment [1], nutrition [2], and resistance to infection [3]. 
In this study, we compared the nasal cavity microbiota 
and the oral cavity microbiota from healthy adult 
humans. Until recently, the bacterial community of the 
healthy human nasal cavity had not been characterized 
by culture-independent methods [4,5]. However, many 
studies have characterized the microbiota of the healthy 
human nares [6-11], which are adjacent and anterior to 
the nasal cavity. Staphylococcus aureus carriage in the 
nares is linked to increased risk of S. aureus infection in 
other body sites [12,13]. Further, antagonism by and 
competition with other members of the nares microbiota 
seem to influence S. aureus nares carriage [12]. Al- 
though adjacent to the nares, the nasal cavity is distinct 
from the nares with a different type of epithelium, a 
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non-keratinized stratified squamous epithelium that 
transitions to a typical respiratory epithelium — ciliated 
pseudostratified columnar epithelial cells and mucus- 
producing goblet cells [14]. In contrast, the nares have a 
keratinized, stratified squamous epithelium with hairs 
and sebaceous glands. Relatively few studies have investi- 
gated the bacterial community composition of the nasal 
cavity in healthy humans. In this study, we sought to ex- 
pand our knowledge of the healthy human nasal cavity 
microbiota and compare the nasal cavity microbiota to 
the oral cavity microbiota in the same subjects. 

Methods 

Subject recruitment and characteristics 

This study was approved by the University of Michigan 
Institutional Review Board. All subjects provided written 
informed consent. Twelve adults patients were recruited 
from a tertiary care otolaryngology clinic (Additional file 1: 
Table SI). Exclusion criteria were patients who had acute 
or chronic sinusitis and patients who were taking antibi- 
otics or oral steroids for any reason. 



© 2014 Bassis et al.; licensee BloMed Central ttd. This is an Open Access article distributed under the terms of the Creative 
Commons Attribution ticense (http://creativecommons.Org/iicenses/by/4.0), which permits unrestricted use, distribution, and 
reproduction In any medium, provided the original work is properly credited. The Creative Commons Public Domain 
Dedication waiver (httpy/creativecommons.org/publlcdomain/zero/1.0/) applies to the data made available In this article, 
unless otherwise stated. 



Bassis ef al. Microbiome 2014, 2:27 
http://www.microbiomejoumal.eom/content/2/1/27 



Page 2 of 5 



Sampling 

4N6 DNA flocked swabs (Cat. No. 3520CA, Copan 
Diagnostics Inc., Murrieta, CA, USA) were used to sam- 
ple all sites. The nasal cavity was sampled by inserting 
the swab into the nasal passage between the septum and 
middle turbinate, taking care to avoid contact to the 
nares. The dorsum of the tongue and buccal mucosa 
were sampled with separate swabs. The samples were 
transferred directly into the Eppendorf tubes provided 
with the swab and stored on ice and then at -20°C until 
DNA isolation. 

DNA isolation 

DNA was isolated from the swabs with a PowerSoil DNA 
isolation kit (Mo Bio Laboratories, Inc., Carlsbad, CA, 
USA) according to the manufacturer's instructions except 
that 2 min of bead beating using the 'Homogenize' setting 
of a Mini-BeadBeater-8 (Biospec Products, Bartlesville, 
OK, USA) was done in place of 10 min of vortexing. 

Primary PCR amplification, pooling, and sequencing 

We based our protocol for amplifying and preparing li- 
braries of the V5V3 region of the 16S rRNA-encoding 
gene on HMP 16S Protocol Version 4.2 (http://www. 
hmpdacc.org/doc/16S_Sequencing_SOP_4.2.2.pdf). Each 
20 |il polymerase chain reaction (PCR) reaction contained 
2 |il AccuPrime PCR Buffer II (Invitrogen, Carlsbad, CA, 
USA), 0.15 1^1 AccuPrime Taq DNA Polymerase High 
Fidelity (Invitrogen), 0.2 primer A (CCATCTCATCC 
CTGCGTGTCTCCGACTCAGXXXXXCCGTCAATTCM 
TTTRAGT), 0.2 ^iM primer B (CCTATCCCCTGTGTG 
CCTTGGCAGTCTCAGCCTACGGGAGGCAGCAG), 
and 1 [A DNA for the oral cavity samples or 15.45 [A DNA 
for the nasal samples. The bold portions of primer A and 
primer B are 926R and 357 F, respectively. The region of 
primer A represented by XXXXX is the 5-10 nucleotide 
barcode sequence. The remainder of primer A and primer 
B are the A adapter sequence and the B adapter sequence, 
respectively, required for emPCR and 454 sequencing. 
The PCR was run for 2 min at 95°C followed by 30 cycles 
of 95°C for 20 s, 50°C for 30 s, and 72°C for 5 min. The 
PCR products were purified with AMPure XP (Agencourt 
Bioscience Corporation, Beckman Coulter, Inc., Beverly, 
MA, USA) according to the manufacturer's instructions 
except 0.6x the amplicon volume (10.8 |il) of beads was 
used rather than 1.2x in order to remove more of the 
small products. The purified PCR products were quanti- 
fied with a Quant-iT PicoGreen dsDNA kit (Invitrogen) 
according to the manufacturer's instructions and com- 
bined into a pool with equal amounts of each amplicon. 
The pool was then purified with AMPure XP (Agencourt 
Bioscience Corporation) according to the manufacturer's 
instructions except the volume of beads was 0.6x the 
pool volume. The pool was quantified with a Library 



Quantification Kit for Roche 454 GS Titanium (Kapa 
Biosystems, Inc., Wilmington, MA, USA). A Junior 
emPCR (454 Life Sciences, Roche, Branford, CT, USA) 
was performed, and 454 sequencing was done on a GS 
Junior (454 Life Sciences) according to manufacturer's 
instructions. 

Sequence processing 

Sequences were processed with mothur v. 1.28.0 according 
to the Schloss SOP of November 27, 2012 [15,16]. In sum- 
mary, the sfif file was input to sffinfo, trim.flows was run 
allowing 1 mismatch in the barcode and 2 mismatches in 
the 926R region of the primer, and sequencing error was 
reduced with shhh.flows. With trim.seqs, barcode and pri- 
mer sequences were removed and all sequences less than 
200 bases or with homopolymers longer than eight nucleo- 
tides were discarded. The sequences were aligned to the 
Silva reference alignment [17,18]. In order to compare se- 
quences over the same region of the alignment, we set the 
end position at 27659 and chose a start position that was 
met by 95% of the sequences. With pre.cluster, sequences 
within two base pairs were merged. Chimeras were identi- 
fied with chimera. uchime [19] and removed. The se- 
quences were classified using a modified form of RDP 
training set version 9 (trainset9_032012.pds.tax and train- 
set9_032012.pds.fasta) [20]. Sequences classified as Chloro- 
plast. Mitochondria, Archaea, Eukaryota, or unknown 
kingdom were removed. 

Sequence analysis 

For our initial analysis, we included subjects only if se- 
quences were obtained for a complete set of samples 
(buccal mucosa, tongue, and nasal cavity). Therefore, in 
order to maximize the number of subjects included in 
the study, we decided to subsample 269 sequences from 
each sample. So, subjects were included in our initial 
analysis only if at least 269 sequences were obtained 
from all three sites (nasal cavity, buccal mucosa, and 
dorsal side of tongue). A distance matrix made with dist. 
seqs was used with the average neighbor algorithm to 
group sequences into operational taxonomic units 
(OTUs) with the cluster command. OTUs defined as 3% 
different were used for further analysis. The make. shared 
command was used to produce a table (shared file) of 
the number of sequence reads assigned to each OTU in 
each sample. The shared file was used to calculate 6 
{1-6 similarity index), a metric of community dissimi- 
larity that takes the relative abundances of both shared 
and non-shared OTUs into account [21]. Principle co- 
ordinates analysis (PCoA) was used to visualize the 6 
distance matrix, and analysis of molecular variance 
(AMOVA) [22] was used to test the statistical signifi- 
cance of the differences between bacterial communities 
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of different groups (i.e., nasal cavity communities versus 
buccal mucosa communities). 

Quality assurance 

To ensure that the source of bacterial sequences was not 
the swab itself or the DNA isolation reagents, PCR was 
performed on DNA isolated from an unused swab. To 
confirm that the PCR reagents were not the source of 
bacterial sequences, PCR of the no-template control was 
performed. Neither of these control PCRs yielded prod- 
ucts visible on a gel, indicating that there was minimal 
contamination from the swab or reagents. 

Initial findings 

Sequences obtained 

After sequence processing, we obtained a total of 35,142 
high-quality bacterial 16S rRNA-encoding gene se- 
quence reads from 35 buccal mucosa, tongue (dorsal 
side), and nasal cavity samples from 12 subjects with a 
mean of 1,115 sequences per sample (Additional file 1: 
Table SI). The mean sequence length after sequence 
processing was 268 bases. The number of sequences ob- 
tained per sample from the buccal mucosa microbiota 
ranged from 715-1,684. The number of sequences ob- 
tained per sample from the tongue microbiota ranged 
from 519-1,597. The number of sequences obtained per 
sample from the nasal cavity microbiota ranged from 1- 



1,595. The sff files from which the processed sequences 
were generated were submitted to the SRA (Bioproject: 
PRJNA248297) except for the sff that yielded only one 
sequence (nasal cavity L) (Additional file 1: Table SI). 

Nasal cavity bacterial communities were distinct from the 
bacterial communities on the tongue and buccal mucosa 

For our initial analysis, we subsampled 269 sequences 
from each sample following sequence processing. Ten 
healthy adult subjects had at least 269 sequences from 
each sampling site following sequence processing (nasal 
cavity, dorsum of the tongue, and buccal mucosa) and 
were included in the analysis (Additional file 1: Table SI). 
As illustrated by PCoA of 6 dissimilarities, the bacter- 
ial communities clustered based on body site, with the 
nasal cavity bacterial communities more similar to 
each other than to the oral cavity bacterial communi- 
ties (Figure 1). Based on 9 dissimilarities, the differ- 
ences between the bacterial communities of the nasal 
cavity and each oral cavity subsite were statistically sig- 
nificant (AMOVA p value <0.001 for each compari- 
son). Additionally, although the bacterial communities 
of the two sites in the oral cavity were more similar to 
each other than to the nasal cavity communities, based 
on 6 dissimilarities, the bacterial communities from 
the tongue and the buccal mucosa were distinct (AMOVA 
/-value <0.001). 
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Figure 1 Principal coordinates analysis of nasal cavity, buccal mucosa and tongue microbiota. Principal coordinates analysis was 
performed on a matrix of 9 distances between all samples. OTUs were defined as 3% different. Axis 1, representing 29.8% of the variance, and 
axis 2, representing 19.3% of the variance, are displayed. 
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Figure 2 Bacterial community composition of the nasal cavity, buccal mucosa, and tongue. Two hundred sixty-nine sequences were 
subsampled from each sample. The sequences were then classified as described in the materials and methods, and the relative abundances of 
sequences classified to major taxonomic groups are displayed. 



Composition of the nasal cavity microbiota 

In our study and other recent studies [4,5], the bacterial 
communities of the nasal cavities were dominated by Acti- 
nobacteria, Firmicutes, and in some cases, Proteobacteria 
(Figure 2). Corynebacteriaceae and Propionibacteriaceae 
were the most prevalent families of Actinobacteria in the 
nasal cavity (Figure 2). Between subjects, the levels of 
Corynebacteriaceae varied from 1.5% to 62.8%, and the 
levels of Propionibacteriaceae varied from 0.4% to 42.4% 
(Figure 2, Additional file 2: Table S2). Actinobacteria were 
present at lower levels in the communities at both oral 
cavity sites, but they were from the families Micrococ- 
caceae and Actinomycetaceae. Corynebacteriaceae com- 
posed over 1% of the community in only one oral cavity 
sample. Propionibacteriaceae was undetectable in most 
oral cavity samples, and when it was detected, it never 
exceeded 1% of the community. 

Another striking difference between the nasal cavity and 
oral cavity communities was the composition of the Firmi- 
cutes. The same Staphylococcaceae OTU was present in 
the nasal cavity samples of all subjects, ranging from 2.2%- 
55.0% of the community, while Staphylococcaceae were not 
detected in the oral cavity communities (Additional file 2: 
Table S2). It was not possible to distinguish Staphylococcus 
aureus and Staphylococcus epidermidis in the region of the 
16S rRNA gene covered by our sequences. Generally, in the 
oral cavity, Streptococcaceae was the most abundant Firmi- 
cutes family, and in several subjects, the most abundant 
family overall. Veillonellaceae was also more abundant in 
the oral cavity. 

The overall levels of Betaproteobacteria were not signifi- 
cantly different between nasal cavity communities and oral 
cavity communities (Kruskal-Wallis test,/? value = 0.4557). 
However, the composition of the Betaproteobacteria dif- 
fered between the nasal cavity and the oral cavity: all three 



sites included Neisseriaceae at levels that were not statisti- 
cally different (Kruskal-Wallis test, p value = 0.4543), but 
the nasal cavity contained higher levels of Comamonada- 
ceae (Kruskal-Wallis test, p value = 0.0002) and Burkhol- 
deriales incertae sedis (Kruskal-Wallis test, p value = 
0.0004). 

Future directions 

With its proximity to the sinuses — the maxillary sinuses 
communicate with the nasal cavity through 1 to 2-mm 
ostia — the nasal cavity microbiota might be a useful 
proxy for the less accessible sinus microbiota. Signatures 
of sinusitis in the sinus microbiota have recently been 
identified [23]. To determine if specific changes in the 
nasal cavity microbiota also accompany sinusitis and 
could be used as an indicator of sinus infections, it will 
be necessary to sample both the nasal cavity and the 
sinus in the same subjects with and without sinusitis. 

Availability of supporting data 

The sff files have been deposited in the SRA (Bioproject: 
PRJNA248297). 

Additional files 



Additional file 1: Table 51. Number of sequences obtained from each 
sample and subject characteristics. 

Additional file 2: Table 52. Relative abundances of bacterial families in 
each sample based on classification of 16S rRNA-encoding gene sequence 
reads. 
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